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Abstract. We report on photometric, spectroscopic and polarimetric monitoring of the optical and near-infrared (NIR) afterglow 
of GRB020405. Ground-based optical observations, performed with 8 different telescopes, started about 1 day after the high- 
energy prompt event and spanned a period of ~ 10 days; the addition of archival HST data extended the coverage up to ~ 150 days 
after the GRB . We report the first detection of the afterglow in NIR bands. The detection of Balmer and oxygen emission lines in 
the optical spectrum of the host galaxy indicates that the GRB is located at redshift z = 0.691. Fe n and Mg n absorption systems 
are detected at z = 0.691 and at z = 0.472 in the afterglow optical spectrum. The latter system is likely caused by absorbing 
clouds in the galaxy complex located ~2" southwest of the GRB020405 host. Hence, for the first time, the galaxy responsible 
for an intervening absorption line system in the spectrum of a GRB afterglow is spectroscopically identified. Optical and NIR 
photometry of the afterglow indicates that, between 1 and 10 days after the GRB, the decay in all bands is consistent with a 
single power law of index a = 1.54+0.06. The late-epoch VLT ./-band and HST optical points lie above the extrapolation of this 
power law, so that a plateau (or "bump") is apparent in the VRIJ light curves at 10-20 days after the GRB. The light curves at 
epochs later than day ~20 after the GRB are consistent with a power-law decay with index a' = 1.85+0. 1 5 . While other authors 
have proposed to reproduce the bump with the template of the supernova (SN) 1998bw, considered the prototypical 'hypernova', 
we suggest that it can also be modeled with a SN having the same temporal profile as the other proposed hypernova SN2002ap, 
but 1.3 mag brighter at peak, and located at the GRB redshift. Alternatively, a shock re-energization may be responsible for 
the rebrightening. A single polarimetric i?-band measurement shows that the afterglow is polarized, with P = 1.5+0.4 % and 
polarization angle 6 = 172°±8°. Broad-band optical-NIR spectral flux distributions show, in the first days after the GRB, a 
change of slope across the J band which we interpret as due to the presence of the electron cooling frequency v c . The analysis 
of the multiwavelength spectrum within the standard fireball model suggests that a population of relativistic electrons with index 
p ~ 2.7 produces the optical-NIR emission via synchrotron radiation in an adiabatically expanding blastwave, with negligible 
host galaxy extinction, and the X-rays via Inverse Compton scattering off lower-frequency afterglow photons. 
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1. Introduction 

GRB020405 was detected by the Third Interplanetary Network 
(IPN) on 2002 April 5.02877 UT with a duration of -40 s and 
localized to an error box of 75 square arcmin size. In the 25- 
100 keV band it had a total fluence of ~3 xl0~ 5 erg cirT 2 and a 
peak flux of ~ 10~ 6 erg crrT 2 s (Hurley et al. 2002). This 
GRB was also observed by the GRBM onboard BeppoSAX, 
with a duration of -60 s in the 40-700 keV band and a 50- 
700 keV fluence of ~4xl0~ 5 erg cirT 2 (Price et al. 2003). 

Approximately 18 hours after the GRB, Price et al. (2002a, 
2003) detected a relatively bright (R ~ 18.5) source within the 
IPN error box which was not present on the DSS-II red plate. 
This Optical Transient (OT), located at coordinates (J2000) RA 
= 13 h 58 m 03-U2; Dec = -31° 22' 22'.'2 (with an error of 073 
along both directions), was confirmed by subsequent observa- 
tions and identified as the afterglow of GRB020405 (Castro- 
Tirado et al. 2002; Palazzi et al. 2002; Hjorth et al. 2002; Price 
et al. 2002b; Gal- Yam et al. 2002a; Covino et al. 2002a,b). 
Optical spectroscopy allowed Masetti et al. (2002a) and Price 
et al. (2003) to determine the redshift of the GRB, z = 0.691 
(see also Sect. 3.3). 

A counterpart to the GRB has been detected also at ra- 
dio and X-ray wavelengths (Mirabal et al. 2002; Berger et al. 
2003). By fitting to the optical light curves a power law F (f) 
oc t~ a , typical of GRB afterglows, Price et al. (2002b) and 
Covino et al. (2002b) found a ~ 1.26 and a = 1.52+0.12, re- 
spectively, from observations performed within 3 days after the 
GRB. By using a data set spanning ~5 days, Price et al. (2003) 
found instead that a better fit to the early optical afterglow light 
curve was obtained with a smoothly broken power law with 
decay indices a\ ~ 0.9 and a% ~ 1.9 before and after a break 
which occurred ~1.7 days after the GRB, respectively. On the 
other hand, Bersier et al. (2003) found no evidence of a break 
between 1.24 and 4.3 days after the GRB in their /?-band data, 
which followed a single power-law decay with index a ~ 1.7. 

Polarimetric measurements by Covino et al. (2003) and 
Bersier et al. (2003) suggest early rapid variability of the per- 
centage of linear polarization between values of ~1.5% and 
-10%. 

Price et al. (2003), using late-time HST observations cov- 
ering the time interval between 20 and 140 days after the high- 
energy prompt event, found a "red bump", or flattening, in the 
optical light curves of the OT and suggested that this could 
be due to an emerging supernova (SN) component, which they 
modeled with the template of SN1998bw located at the redshift 
of GRB020405 and further dimmed by about 0.5 mag. A simi- 
lar assumption was made by Dado et al. (2002) in their model- 
ing of the GRB020405 optical afterglow within the cannonball 
picture. 

Optical observations also revealed the presence of a rela- 
tively large (about 2"xl", elongated in the N-S direction) neb- 

Send offprint requests to: N. Masetti, masetti@bo . iasf . cnr . it 

* Based on observations made with ESO telescopes at Paranal and 
La Silla Observatories under programme ID 69.D-0701, with the tele- 
scopes TNG, WHT and JKT, operating in the Spanish Observatorio 
del Roque de los Muchachos of the Instituto de Astroffsica de 
Canarias, and with the lm telescope of SO in NainiTal, India 



ulosity located ~2" southwest of the OT and interpreted as its 
putative host galaxy (Hjorth et al. 2002). 

In this paper we report on optical imaging, spectroscopy 
and polarimetry, along with the first detection in the near- 
infrared (NIR) bands of the GRB020405 afterglow. The data 
were acquired, in the framework of the GRACE 1 collaboration, 
at the ESO telescopes of Cerro Paranal and La Silla (Chile) 
starting on 6 April 2002, i.e., ~1 day after the high-energy 
prompt event. Observations from the 1 -metre SO telescope, 
WHT and TNG are also included, along with the reanalysis 
of the public HST observations of this GRB. 

The paper is organized as follows: Sect. 2 describes the op- 
tical and NIR observations and the data analysis; the results 
are reported in Sect. 3 and discussed in Sect. 4; in Sect. 5 we 
report our conclusions. Throughout the paper we assume a cos- 
mology with H = 65 km s _1 Mpc~', Q A = 0.7 and Q m = 0.3; 
also, when not otherwise stated, uncertainties will be reported 
at lcr confidence level, and upper limits at 3<x confidence level. 

2. Observations and data reduction 

2.1. Optical and NIR photometry 

Optical BVRI observations in the Johnson-Cousins photomet- 
ric system were performed with YLT-Melipal plus FORS1, 
VLT-Yepun plus FORS2, NTT plus SUSI2 and EMMI, and 
1.54-metre Danish telescope plus DFOSC. 

FORS1 is equipped with a 2048x2048 pixels Tektronix 
CCD which covers a 6f8x6!8 field in the standard resolution 
imaging mode with a scale of 072 pix~'; FORS2 is equipped 
with a mosaic of two 2048x4096 pixels MIT CCDs, usually 
operating in 2x2 standard binning mode, thus covering a field 
with size 6.'8x6.'8 with a spatial resolution of 0725 pix -1 in the 
standard resolution imaging mode. SUSI2 incorporates a mo- 
saic of two 2048x4096 pixels EEV CCDs which cover a field 
of 5:5x5:5, corresponding to 0708 pix -1 ; EMMI (Red Arm) 
carried until 14 May 2002 a 2048x2047 pixels Tektronix CCD 
with a scale of 0727 pix~' and a field coverage of 9'15x8'6 in 
size. DFOSC is equipped with a 2048x2048 pixels CCD cover- 
ing a 13:7x1377 field, thus securing a spatial resolution of 0739 
pix -1 . 

Optical observations were also obtained on 5 and 6 April 
2002 (/ band) at the 1.0m Sampurnanand telescope of SO lo- 
cated in Nainital (India), and on 6 April 2002 at the Canary 
Islands (Spain) with the 4.2m WHT plus PFIP (UBVRI bands) 
and with the 3.58m TNG plus DOLoReS (V band). The 1.0m 
SO telescope carries a 2048x2048 pixels CCD, with a 13'xl3' 
field of view; in its standard 2x2 binning mode it has a spa- 
tial resolution of 0776 pix~'. The imaging camera PFIP car- 
ries two 2100x4200 EEV CCDs which cover a field of view 
of 16:2x16:2, giving a plate scale of 0724 pix~'; the spec- 
trophotometer DOLoReS carries a 2048x2048 pixels Loral 
CCD which can image a field of 9:5x9:5 with a scale of 07275 
pix -1 . 

NIR imaging was obtained at ESO NTT with SofI (/, H and 
K s bands) and at WW-Antu with ISAAC (J s band). The SofI in- 

1 GRB Afterglow Collaboration at ESO: see the web page 
http : / /zon . wins . uva . nl/grb/ " grace 
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frared spectrograph and imaging camera works in the 0.9-2.5 
/im NIR range by using a Hawaii 1024x1024 pixel HgCdTe 
array. In the small-field imaging mode the plate scale is 07144 
pix -1 and the corresponding field of view is 2'Ax2'A. ISAAC 
is equipped, in the 0.9-2.5 fim range, with a Rockwell Hawaii 
1024x1024 pixel HgCdTe array which has a scale of a.' 148 
pix 1 and secures imaging with a field of view of 2'5x2'5. The 
J s filter is centered at 1 .24 fim and has a full width at half max- 
imum of 0.16 fim; its overall response is 25% lower than the 
J filter. This is confirmed by comparing the nightly zero-point 
coefficients of / and J s frames in our data set. We thus reported 
the ISAAC J s magnitudes to the standard / band by increasing 
the observed 7 s -band flux densities by 25%. In order to allow 
for sky subtraction, the total integration time of each NIR point- 
ing was split into dithered images of 15 s each when using Sofl, 
and of 30 s each during ISAAC observations. In both cases the 
dithering was 40" between consecutive images. 

The complete log of our optical and NIR imaging observa- 
tions is reported in Table 1 . 

Optical images were bias-subtracted and flat-fielded with 
the standard cleaning procedure. In some cases (especially at 
late epochs), frames taken on the same night in the same band 
were summed together in order to increase the signal-to-noise 
ratio. In Fig. 1 we report the R-band image of the field of the 
GRB020405 counterpart obtained on 6 April 2002 with the 
1.54-metre Danish telescope plus DFOSC. 

Since the close environment of the afterglow is quite 
crowded (see Figs. 3 and 6), we chose standard Point Spread 
Function (PSF) fitting, rather than simple aperture, photome- 
try. We used the DAOPHOT II image data analysis package 
PSF-fitting algorithm (Stetson 1987) running within MIDAS 2 . 
A two-dimensional Gaussian profile with two free parameters 
(the half width at half maxima along x and y coordinates of 
each frame) was modeled on at least 5 unsaturated bright stars 
in each image. The errors associated with the measurements 
reported in Table 1 represent statistical uncertainties obtained 
with the standard PSF-fitting procedure. 

The BVRI zero-point calibration was performed by using 
several standard fields (Landolt 1992) taken under photometric 
conditions at several ESO telescopes. The single Gunn i image 
(6 April 2002) was calibrated using the /-band secondary stan- 
dards, given that the widths, the reference wavelengths and the 
flux density normalizations of the two filters are very similar 
(Fukugita et al. 1995). However, in order to account for small 
differences between the filters, we added in quadrature a 3% 
uncertainty to the /-band magnitude error obtained from the 
Gunn i observation. 

Concerning the single f/-band measurement taken at WHT, 
no standard field was acquired on the same night due to non- 
photometric conditions of the sky. The calibration in this band 
was then performed on 10 May 2002 at the 1.0-metre JKT, lo- 
cated in the Canary Islands and equipped with a SITe2 CCD 
which has 2048x2048 pixels and an image scale of 0.33 pix~'. 
Observations of the GRB field and of some Landolt (1992) 

2 MIDAS (Munich Image Data Analysis System) is developed, dis- 
tributed and maintained by ESO (European Southern Observatory) 
and is available at http: //www. eso . org/projects/esomidas/ 



fields were acquired in U and B to determine the CCD color 
term. 

We then selected 6 stars of different brightness in the 
GRB020405 field and we used them to determine the OT op- 
tical BVRI magnitudes (see also Simoncelli et al. 2002). Two 
additional stars were used for the calibration in the U band. All 
these stars are indicated in Fig. 1 and their magnitudes are listed 
in Table 2. We find this magnitude calibration to be accurate to 
within 3%. We note that the B and R magnitudes of field stars 
as reported in the USNO-A2.0 catalog 3 differ by ~0.6 mag with 
respect to our calibration, which is however consistent with that 
of Bersier et al. (2003) and Price et al. (2003). 

Public HST/WFPC2 observations of GRB020405 taken in 
three filters (F555W, F702W and F814W) have been included 
in our analysis (see also Price et al. 2003). All HST data were 
retrieved from the multimission archive 4 after "On-The-Fly" 
calibration had been performed. Dithered images were then 
combined and cosmic -ray cleaned using the drizzle method 
(Fruchter & Hook 2002). The images were drizzled onto a final 
output grid with pixels half the size of the native WPFC2 pixels 
using a pixfrac parameter of 0.7. 

The images acquired with filters F555W and F702W on 
August 23, 2002, i.e., 140 days after the GRB, do not show 
substantial pointlike emission at the OT position; therefore they 
were used to evaluate the host galaxy contribution underlying 
the afterglow. This could not be done for the F814W data as 
no images were acquired with this filter during the last HST 
visit. Therefore, in order to estimate the magnitude of the host 
underlying the OT position in the F814W frames, we created a 
PSF using the Tiny Tim software (Krist 1993). This PSF was 
subsampled and convolved to mimic the effect of the drizzle. 
The PSF was then subtracted from the position of the OT in 
the June 2002 F814W HST images to the level where the re- 
maining surface brightness was entirely due to the host. These 
PSF-subtracted images were then used as a template to deter- 
mine and subtract the host galaxy contribution in all previous 
(space- and ground-based) observations made in that band. The 
accuracy of this method was checked by applying it to the June 
2002 images acquired with filters F555W and F702W, and by 
comparing the results with those obtained using the August 
2002 data. We found that the host magnitudes measured at the 
OT position with the two methods (i.e., direct host contribution 
measurement on one hand, and PSF modeling plus subtraction 
on the other) coincide within the uncertainties. This made us 
confident that the host contribution, as derived with the PSF- 
subtraction method, was correct and could be applied to the 
F814W images. 

OT magnitudes were then measured on each available HST 
image with standard aperture photometry, and were calibrated 
to the Vega system using the published zeropoints for the 
WFPC2 (Biretta et al. 2002). We assumed that the conversion 
from the WFPC2 to the Johnson-Cousins system did not induce 
significant errors on the photometry, as also stated in Fukugita 
et al. (1995). Indeed, the resulting OT magnitudes, reported in 
Table 3, are consistent with the flux values independently re- 

3 available at http: //www. nofs .navy .mil/ 

4 see http: //archive. stsci.edu/ 
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Table 1. Journal of the GRB020405 afterglow ground-based optical and NIR photometric observations. Magnitudes are not 
corrected for Galactic interstellar absorption 



Mid-exposure 


Telescope 


Instrument 


Filter 


Total exposure 


Seeing 


Magnitude 


time (UT) 








time (s) 


(arcsec) 




2002 Apr 5.812 


lm SO 


CCD 


/ 


2x600 


2.6 


19.30±0.06 


2002 Apr 6.069 


TNG 


DOLoRes 


V 


2x300 


1.8 


20.66+0.02* 


2002 Apr 6.013 


WHT 


PFIP 


R 


300 


2.3 


20.17+0.03 


2002 Apr 6.109 


WHT 


PFIP 


U 


900 


2.2 


20.48+0.06 


2002 Apr 6.118 


WHT 


PFIP 


B 


600 


2.5 


21.21+0.03 


2002 Apr 6.124 


WHT 


PFIP 


V 


300 


2.0 


20.73+0.03 


2002 Apr 6.129 


WHT 


PFIP 


R 


300 


2.8 


20.28+0.04 


2002 Apr 6.134 


WHT 


PFIP 


I 


300 


3.2 


19.61+0.08 


2002 Apr 6.136 


1.54D 


DFOSC 


R 


600 


1.4 


20.38+0.02 


2002 Apr 6.186 


Melipal 


FORS1 


R 


60 


0.7 


20.48+0.01 


2002 Apr 6.197 


Melipal 


FORS1 


R 


30 


0.8 


20.46+0.0 1 


2002 Apr 6.211 


1.54D 


DFOSC 


R 


600 


1.4 


20.44+0.02 


2002 Apr 6.220 


1.54D 


DFOSC 


V 


600 


1.4 


20.81 ±0.02 


2002 Apr 6.220 


NTT 


SofI 


H 


300 


1.0 


18.07+0.07 


2002 Apr 6.224 


NTT 


Sofl 


J 


300 


1.0 


18.69+0.06 


2002 Apr 6.230 


1.54D 


DFOSC 


B 


600 


1.4 


21.31+0.03 


2002 Apr 6.230 


NTT 


Sofl 


K s 


600 


0.9 


17.34+0.06 


2002 Apr 6.234 


Melipal 


FORS1 


R 


20 


0.8 


20.53+0.02 


2002 Apr 6.246 


1.54D 


DFOSC 


Gunn i 


600 


1.1 


19.88+0.02 


2002 Apr 6.797 


lm SO 


CCD 


I 


5x600 


2.5 


>20.0 


2002 Apr 7.132 


1.54D 


DFOSC 


R 


600 


1.2 


21.20+0.04 


2002 Apr 7.240 


1.54D 


DFOSC 


R 


600 


0.9 


21.30+0.06 


2002 Apr 7.249 


1.54D 


DFOSC 


V 


600 


1.1 


21.79+0.05 


2002 Apr 7.266 


1.54D 


DFOSC 


B 


600 


1.2 


22.29+0.07 


2002 Apr 7.270 


NTT 


Sofl 


H 


600 


0.8 


18.90+0.07 


2002 Apr 7.278 


NTT 


Sofl 


J 


600 


0.8 


19.55+0.06 


2002 Apr 7.281 


Melipal 


FORS1 


V 


2x30 


0.7 


21.95+0.04 


2002 Apr 7.294 


NTT 


Sofl 


K s 


1800 


0.7 


18.18+0.06 


2002 Apr 7.381 


1.54D 


DFOSC 


R 


600 


1.2 


21. 32+0.06 


2002 Apr 8.280 


NTT 


SUSI2 


V 


10x60 


0.8 


22.49+0.02 


2002 Apr 8.294 


NTT 


SUSI2 


R 


10x60 


0.8 


22.09+0.03 


2002 Apr 8.310 


NTT 


SUSI2 


I 


5x60+3x300 


0.8 


21.44+0.02 


2002 Apr 9.205 


NTT 


EMMI 


R 


2x300 


1.0 


22.23+0.04 


2002 Apr 9.194 


NTT 


EMMI 


V 


2x300 


1.0 


22.76+0.06 


2002 Apr 9.215 


NTT 


EMMI 


I 


2x300 


1.0 


21.70+0.06 


2002 Apr 12.097 


Yepun 


FORS2 


V 


3x180 


0.7 


23.59+0.04 


2002 Apr 12.107 


Yepun 


FORS2 


R 


3x180 


0.7 


23.09+0.04 


2002 Apr 15.238 


Melipal 


FORS1 


V 


2x180 


0.7 


23.99+0.07 


2002 Apr 15.244 


Melipal 


FORS1 


R 


2x180 


0.6 


23.36+0.04 


2002 Apr 15.250 


Melipal 


FORS1 


1 


2x180 


0.6 


22.53+0.06 


2002 Apr 27.193 


Antu 


ISAAC 


Js 


3600 


0.5 


22.4+0.1 


2002 May 8.193 


Antu 


ISAAC 


Js 


3600 


0.6 


>22.7 



*A preliminary analysis of this observation was presented in Palazzi et al. (2002) 



ported by Price et al. (2003). The only discrepancy is found in 
the F8 14W filter data and arises from the fact that those authors 
assumed no OT contribution in the June 2002 HST image ac- 
quired with that filter, so that was used directly as the template 
for the host galaxy subtraction. In fact, if we subtract, for that 
HST filter, the June 2002 magnitude of the OT plus underlying 
host from those obtained in the previous visits, we obtain flux 
density values consistent with those reported in Table 3 of Price 
et al. (2003). 



The reduction of the NIR images was performed with IRAF 
and the STSDAS packages 5 . Each image was reduced by first 



5 IRAF is the Image Analysis and Reduction Facility made 
available to the astronomical community by the National Optical 
Astronomy Observatories, which are operated by AURA, Inc., under 
contract with the U.S. National Science Foundation. STSDAS is dis- 
tributed by the Space Telescope Science Institute, which is operated by 
the Association of Universities for Research in Astronomy (AURA), 
Inc., under NASA contract NAS 5-26555. 
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Table 2. Optical magnitudes of selected GRB020405 field stars as indicated in Fig. 1. Values are not reported for cases in which 
the stars are either saturated or too faint for reliable calibration 



Star 


U 


B 


V 


R 


/ 


1 




1Q 01 

1 ./ . 1 -!—\J.\J ± 


± (J . \J S -!—\J ,\J ± 


18 17+0 01 


17 ^5+0 01 


2 




19.84±0.01 


19.01±0.01 


18.48±0.01 


17.87±0.01 


3 




19.34±0.01 


18.85±0.01 


18.50±0.02 


18.12+0.01 


4 




20.09±0.02 


19.49±0.01 


19.14±0.02 


18.72±0.01 


5 




18.18+0.01 


17.92±0.01 


17.72±0.01 


17.45±0.01 


6 




20.45±0.02 


19.35±0.01 


18.67±0.01 


18.00±0.01 


A 


15.76+0.01 


16.00±0.02 








B 


18.35±0.03 


17.74±0.02 
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Fig. 1. /?-band image of the field of 
GRB020405 acquired with the 1.54- 
metre Danish telescope plus DFOSC 
on 2002 April 6.136 UT at ESO-La 
Silla. The OT is at the centre of the 
image, indicated with the tick marks. 
North is at top, East is to the left; the 
field size is about 3'x3'. Numbers and 
letters indicate the reference stars with 
magnitudes reported in Table 2 



subtracting a mean sky, obtained from the median of a num- 
ber of images acquired just before and after each processed 
frame. Before the frames were used for sky subtraction, stars in 
them were eliminated by a background interpolation algorithm 
(imedit) combined with an automatic "star finder" (daofind). 
Then, a differential dome flatfield correction was applied; next, 
the telescope dithering was measured from the offsets of field 
objects in each image and the images were averaged together 
using inter-pixel shifts. Also for the NIR images we used PSF- 



fitting photometry following the same procedure adopted in the 
optical. We calibrated the NIR photometry using stars selected 
from the NICMOS Standards List (Persson et al. 1998) and de- 
termined the NIR zero-point coefficients for each observation. 
The standard stars were observed in five positions on the de- 
tector, and their images were reduced in the same way as those 
of the GRB field. Formal photometric accuracy based only on 
the standard stars observations is typically better than 3%. Our 
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calibration is fully consistent with the results of the 2MASS 
survey 6 (Skrutskie et al. 1997) for the GRB020405 field. 

Next we evaluated the Galactic absorption in the optical 
and NIR bands along the direction of GRB020405 using the 
Galactic dust infrared maps by Schlegel et al. (1998); from 
these data we obtained a color excess E(B — V) = 0.055. By 
applying the relation by Cardelli et al. (1989), we derived Au 
— 0.27, Ab — 0.22, Ay — 0.17, Ar — 0.14, A 7 = 0.10, Aj = 0.05, 
A H = 0.03 and A Ks = 0.02. 

2.2. Optical spectroscopy 

Two series of spectra with total exposure times of 30 minutes 
and 60 minutes were obtained on 2002 April 6.213 and 7.310 
UT, respectively, with VLT-Melipal+FORSl. The first series 
was acquired using Grism 1501 plus order separator OG590, 
which avoids overlapping of spectral orders over a given wave- 
length; this limited the spectral range to 6000-9000 A. For the 
second series the Grism 300V with no order separator was 
used: this allowed a wider spectral coverage (3500-9000 A) and 
a spectral resolution higher by about a factor of two. The slit 
width was 1" for both observations and these setups secured 
final dispersions of 5.5 A pix~' for the first series and of 2.6 
A pix~' for the second. In the second spectroscopic pointing 
the slit was rotated by about 40° with respect to the N-S direc- 
tion in order to include both the OT and the galaxy located ~2" 
southwest of it (Hjorth et al. 2002). 

The spectra, after correction for flat-field and bias, were 
background subtracted and optimally extracted (Home 1986) 
using IRAF. He-Ne-Ar and Hg-Cd lamps were used for wave- 
length calibration; both spectroscopic runs were then flux- 
calibrated by using the spectroscopic standards LTT 6248 and 
Hiltner 600 (Hamuy et al. 1992, 1994) for the April 6 and 7 ob- 
servations, respectively. Finally, spectra taken within the same 
night were stacked together to increase the S/N ratio. The cor- 
rectness of the wavelength and flux calibrations was checked 
against the position of night sky lines and the photometric data 
collected around the epoch at which the spectra were acquired, 
respectively. The typical errors were 0.3 A for the wavelength 
calibration and 10% for the flux calibration. 

2.3. Optical polarimetry 

Linear polarimetry in the R band was accomplished on 2002 
April 6.236 UT at VLT-Me/i/?a/+FORS 1 under an average see- 
ing of 0'.'8. Data were obtained by using a Wollaston prism and 
a half-wavelength, phase-retarder plate. The Wollaston prism 
separates the incident light into an ordinary and an extraor- 
dinary component, while the phase-retarder plate determines 
which of the Stokes parameters is measured (U or Q). For each 
image, a mask producing 22"-wide parallel strips was used to 
avoid overlap of the ordinary and extraordinary components. 
The observation therefore consisted of four exposures centered 
at the position of the OT, with the phase-retarder plate at 0°, 
22 °5, 45° and 67 °.5. Each angle was imaged with an expo- 
sure time of 300 s, for a total exposure time of 1200 s over the 

6 available at http : //www. ipac . caltech. edu/2mass/ 



Table 3. HST-WFPC2 optical magnitudes of the GRB020405 
OT. The values are corrected for the host galaxy contribution, 
but not for the Galactic absorption 



Mid-exposure date (UT) 


Filter 


Magnitude 


2002 Apr 24.272 


F555W 


25.05±0.04 


2002 May 5.631 


F555W 


25.83±0.05 


2002 Jun 2.711 


F555W 


27.25±0.15 


2002 Aug 23.246 


F555W 


>27.5 


2002 Apr 28.452 


F702W 


93 QQ+0 02 


2002 May 1.620 


F702W 


24.31 ±0.02 


2002 May 3.625 


F702W 


24.38±0.02 


2002 Jun 1.614 


F702W 


25.91±0.12 


2002 Aug 23.451 


F702W 


>26.7 


2002 Apr 26.275 


F814W 


23.13+0.01 


2002 May 1.486 


F814W 


23.45±0.02 


2002 Jun 9.560 


F814W 


25.42±0.15 



four angles. Image reduction and analysis was performed as de- 
scribed in Sect. 2.1 for the optical photometry. No significant 
variation among the values of the Q and U parameters of indi- 
vidual field stars was noted, therefore we computed the average 
values of these quantities and subtracted them from the corre- 
sponding parameters of the OT to remove the instrumental and 
(local) interstellar polarization. We also checked that the po- 
larization of the selected field stars did not systematically vary 
with their position on the CCD or with their magnitude. The 
polarization angle 6 was calibrated by using the polarized stan- 
dard Velal_95 (Szeifert 2001). The unpolarized standard WD 
1620-315 (Szeifert 2001) was used to check the instrumental 
polarization, which was found to be negligible. 

3. Results 

The optical (see Fig. 1) and NIR images show a point-like 
source at the OT position indicated by Price et al. (2002a, 
2003). The transient behaviour observed in both spectral ranges 
(see Table 1) confirms the afterglow nature of the source. 

3.1. Light curves 

In the analysis of the optical and NIR light curves of the 
GRB020405 we added to our data set the points published by 
Covino et al. (2003), whose photometry can be connected with 
ours. We first considered the ground-based data. These may be 
contaminated by the underlying host galaxy complex. Thus, in 
order to correct for its contribution to the VRI magnitudes of 
the afterglow, we followed two approaches: (i) we fit the light 
curves with a single power law plus constant (F oc r" + const); 
(ii) we measured, on the latest available (OT-subtracted) HST 
images, the contribution of the host galaxy complex within 
an aperture radius matching the ground-based telescopes PSFs 
and subtracted it from our photometric measurements. The two 
methods gave results which are consistent within the errors, and 
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Fig. 2. VRIJ light curves of 
the GRB020405 afterglow. 
Different symbol styles indi- 
cate different bands. Data are 
corrected for the underlying 
host galaxy contribution and for 
Galactic absorption. The light 
curves were rescaled in flux for 
clarity using the factors indi- 
cated in the figure. A flattening 
between day 10 and 20 after the 
GRB, more evident in the RIJ 
bands, is detected. The light 
curves are best fitted using a 
power law with a = 1.54 (up to 
10 days after the GRB) and a 
power law with a' = 1.85 (after 
day 20) 



showed that the host contribution to the OT luminosity was of 
the order of few percent (and thus negligible if compared with 
the uncertainties) in the first days of observation. 

Concerning the /-band data, in the VLT-ISAAC image of 8 
May 2002 we did not detect a point-like source at the position 
of the afterglow, despite the depth and good resolution of the 
image. Therefore we evaluated the host galaxy contribution in 
the / band by using aperture photometry centered at the OT 
position and with radius equal to the FWHM of the image PSF. 
This can be conservatively considered as an upper limit to the 
OT flux. 

No late-time HST or ground-based data are available for 
other optical or NIR filters. Therefore, no host-subtraction 
could be done for bands different than VRIJ. Since the host 
contribution could be evaluated only for these four bands, only 
the corresponding light curves are reported in Fig. 2. 

The host-subtracted VRI light curves between 1 and 10 
days after the GRB decay according to single power laws 
whose temporal indices are not significantly different from fil- 
ter to filter. Their average value is a = 1.54+0.06. The B light 
curve, which consists of few, early epoch points, is fitted by the 
same temporal power law. By modeling the light curves using 
a smoothly broken power law as described by Beuermann et 
al. (1999) we obtain a worse fit (i.e., with a larger reduced^ 2 ), 
indicating that no temporal break is apparent in our optical data 
set. 

While the /-band light curve is reasonably fitted with a sin- 
gle power law with index compatible with that measured for the 
optical data, the H and K s data appear to decay with a shallower 
slope, a ~ 1.3. Considering only the NIR data acquired in the 
first two days (which are presumably least contaminated by the 



host galaxy), we get an average decay slope ar NIR = 1.28+0.03, 
substantially shallower than the one measured in the optical. 

As already noted for this (Price et al. 2003) and for other, 
well sampled, afterglows (e.g., GRB010222: Masetti et al. 
2001), a single or smoothly broken power-law fit does not for- 
mally give an acceptable fit to the VRI data. However, if we 
quadratically add a constant uncertainty of less than 3% to all 
VRI points in our data set, the single power-law model gives 
an acceptable fit (with^- 2 , ~ 1) for the same a obtained before. 
This formal discrepancy between the data and the model might 
be due to instrumental differences (e.g., the use of different de- 
tectors and different telescopes) or to intrinsic short time scale 
variations of the afterglow light curve (see e.g., Masetti et al. 
2000; Holland et al. 2002; Jakobsson et al. 2002; Holland et 
al. 2003), possibly due to inhomogeneities of the circumburst 
medium (Wang & Loeb 2000). 

Between day 10 and 20 after the GRB, the OT decay flat- 
tens in all bands but more evidently in the RIJ bands, thus 
producing the red bump first noted by Price et al. (2003); 
thereafter, the flux decreases at a faster rate (Fig. 2). We fit- 
ted the data after day 20 to a single power law with index 
nearly independent of the wavelength and of average value a' = 
1.85+0.15. Interestingly, the 8.45 GHz light curve of the radio 
afterglow (Berger et al. 2003) closely resembles this late-time 
optical behaviour. 

Finally we note that the HST images in the filters F702W 
and F814W reveal the presence of a transient source located 
in the southern outskirts of galaxy '1' (as defined in Sect. 3.3 
and Fig. 6), at coordinates (J2000) RA = 13 h 58 m 03 S D1; Dec 
= -31° 22' 2477 (J2000), with a conservative uncertainty of 
071 on both. In Fig. 3 the position of this second transient is 
shown. This object faded between the end of April and the be- 
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Fig. 3. HST-WFPC2 images of the OT field acquired with the F702W filter on April 28 (left panel) and August 23 (central 
panel), along with their difference (right panel). Besides the OT, located at the centre of the image, a second transient object (also 
indicated with tick marks in the April 28 image) is apparent ~3" southwest of the OT itself. The images are about 10"xl0" in 
size; North is at top, East is to the left 



ginning of June 2002, and is not detected in the August 2002 
HST visit. It is also not detected in any of the F555W HST im- 
ages. A more detailed description of this source will be given 
in a forthcoming paper. 

3.2. Optical-NIR spectral flux distribution 

We have constructed optical-NIR broad-band spectra at 8 
epochs spanning from ~1 to ~60 days after the GRB and for 
which we have sufficient optical and/or NIR coverage (Fig. 4). 
The data points were corrected for the Galactic absorption and 
converted into flux densities using the tables by Fukugita et al. 
(1995) for the optical and by Bersanelli et al. (1991) for the 
NIR. The spectra of the first two epochs were not corrected for 
the host contribution because this is known for the VRIJ bands 
only; however, given this was quite modest at those epochs (see 
Sect. 3.1), we simply added a 5% error in quadrature to the un- 
certainties on the optical-NIR flux densities. The subsequent 
spectra were also corrected for the host galaxy contribution; 
the errors on the points include the uncertainty on the host pho- 
tometry. 

For the early epoch (within ~10 days after the GRB) broad- 
band spectra, in the cases in which the data in different bands 
were not simultaneous, the flux density was interpolated to the 
reference epoch assuming the best-fit power-law decay (see 
Sect. 3.1). In the last three broad-band spectra the interpola- 
tions to a common epoch were computed using the power-law 
decay slope (a' = 1.85) determined from observations acquired 
after day 20 from the GRB. 

Single power law (F v oc v~P) fits of the first two broad- 
band spectra are unacceptable due to a significant deviation of 
the J-band points. Since we can firmly exclude instrumental 
effects or systematic errors in the /-band calibration, we fitted 
the spectra with broken power laws. The break is found to be 
located at (2.5+0.6)xl0 14 Hz, i.e., in the J band, and the slopes 
are /?nir = 0.65+0.2 and /3 opt = 1.3+0.2, with no significant 
variability of these parameters between the two epochs. The 



optical slope agrees well, within the uncertainties, with those 
determined by Price et al. (2003) and Bersier et al. (2003). 

The broad-band optical spectra of April 8.3, 9.2 and 15.2 
UT were well fitted with single power laws, with spectral 
slopes 1.0+0.2, 0.9+0.3 and 1.35+0.5, respectively; these are 
consistent with the optical slope measured in the first two 
epochs. The OT emission during the "red bump", measured 23, 
33, and 58 days after the GRB, is fitted by a power law of aver- 
age spectral index j3 = 3.5+0.5, in agreement with the findings 
of Price et al. (2003) and substantially steeper than at the ear- 
lier epochs. A remarkable NIR-to-optical spectral curvature is 
observed. 

The optical and NIR colors of the OT of GRB020405 dur- 
ing the first 10 days of monitoring after the high-energy event 
fall in the loci populated by GRB afterglows in the color-color 
diagrams as illustrated by Simon et al. (2001) and by Gorosabel 
et al. (2002). 

3.3. Optical spectra 

The continuum slopes measured for the VLT spectra of April 6 
and 7 are consistent with the results of our broad-band spectral 
analysis reported in the previous Subsection. 

In the OT spectra of April 6 and 7 we detected a number 
of emission lines which we identified with the forbidden tran- 
sitions of [O n] ,0727, [O in] ,U4959,5007, and with the Hp 
line at an average redshift of 0.6908+0.0017. In the spectrum 
of April 7 (Fig. 5) we also possibly detect H r in emission. Our 
redshift measurement is confirmed by the detection of several 
absorption features blueward of 5000 A in the April 7 spec- 
trum. These absorptions are due to rest-frame ultraviolet (UV) 
metallic lines usually present in OT spectra and, more gen- 
erally, in the spectra of high-z objects (see, e.g., Savaglio et 
al. 2003 and references therein). In this wavelength range we 
also detect an absorption system at a lower average redshift z = 
0.4720+0.0011. 

Our line fitting, performed with the SPLOT task within 
IRAF, assumes a Gaussian profile for both emission and ab- 
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Fig. 5. Spectrum of the OT of GRB020405 obtained with VLT-Melipal plus FORS 1 on 7 April 2002. The spectrum was smoothed 
with a Gaussian filter with cr = 3 A (i.e., comparable with the spectral resolution) and corrected for foreground Galactic absorption 
assuming E(B - V) = 0.055. For clarity, we split it into two panels. Numbers mark the positions of the identified absorption lines 
as listed in Table 4, while the question mark indicates an unidentified line. Telluric absorption features at 6870 A and 7600 A are 
indicated with the symbol © 



sorption lines. A conservative error of 6 A and 3 A, comparable 
with the spectral resolution, is associated with line wavelength 
measurements in the spectra of April 6 and 7, respectively. In 
Table 4 are reported, for each line detected in the two spec- 
troscopic runs, the identification along with the observed and 
rest-frame wavelengths and the associated redshift. No signif- 
icant difference in the emission and absorption line redshifts 
is found in either observation for the features at z = 0.691. In 
Table 4 we also list the equivalent widths (EWs) of the iden- 
tified lines computed in the absorber rest frame, i.e., dividing 
the measured value by the factor (1+z). The errors on the EWs 
are computed by assuming for the spectral continuum in prox- 
imity of each line the values corresponding to its 1-cr lower 
and upper bounds. We did not find any significant variation in 
the emission line fluxes between the two nights, as expected 
if the line production sites are gas-rich and star-forming re- 
gions of the GRB host galaxy. The EWs of the only absorp- 
tion lines detected in both nights (those of Ca n) do not vary 
significantly. We also detected a broad absorption feature (the 
observer's frame EW is 5.6±0.5 A) at 4691 A which cannot be 
identified with any known line among those typically detected 



in OT spectra at either of the two detected redshifts. This fea- 
ture is indicated with a question mark in Fig. 5. 

The slit orientation used for the spectroscopic observations 
acquired on April 7 also allowed us to acquire spectra of the 
nearby galaxy first noted by Hjorth et al. (2002). Fig. 6 reports 
a sketch of the setup with which the VET spectrum of April 
7 was obtained: the l"-wide slit was rotated by 40° towards 
East with respect to the N-S direction in order to acquire also 
the spectrum of the galaxy at 2" southwest of the OT, marked 
as '1' in Fig. 6 (see also Hjorth et al. 2002). The slit also in- 
cluded a galaxy (indicated as '2' in the same Figure) located 
~6" southwest of the OT. The spectra of these two galaxies are 
reported in Fig. 7. 

From the detection of [O n] Xil21 and [O m] /15007 emis- 
sion lines in their spectra, galaxies T and '2' appear to be 
both at redshift z = 0.472+0.001, substantially lower than that 
of the GRB. This is the same redshift of one of the absorp- 
tion systems detected in the OT spectrum. Therefore, galaxy 
'1' is most probably responsible for that absorption system. 
Moreover, galaxy ' 1 ' cannot be interacting with the GRB host. 
This result overrides that reported in Masetti et al. (2002b): the 
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Table 4. List of emission and absorption lines identified in the VLT optical spectra of the GRB020405 OT acquired on April 6 
and 7. The numbers in the leftmost column refer to the identifications shown in Fig. 5. The error on line positions is assumed to 
be +6 A for the spectrum of April 6 and ±3 A for that of April 7, i.e., comparable with the corresponding spectral resolution (see 
text). For absorption lines we report the rest-frame EW, and for emission lines their intensity 



Line 


Observed 


Rest frame 


Transition Redshift 


EW res , 


Flux 


number 


wavelength (A) 


wavelength (A) 






(A) 


(xlO -17 erg cirr 2 s~') 


6 April 2002 




6302 


3727 


[On] 


0.6908+0.0016 




9.5+1.5 




6647 


3935 


CanK 


0.6893±0.0015 


1.2+0.3 






6710 


3969 


CanH 


0.6905±0.0015 


1.2+0.3 






8227 


4861 




0.6923+0.0012 




4.0+0.5 




8386 


4959 


[O m] 


0.6910+0.0012 




4.0+0.5 




8472 


5007 


[O m] 


0.6920+0.0012 




11.5+1.5 






Weighted mean 




0.6910+0.0014 






7 April 2002 


1 


3965 


2344 


Fe ii 


0.6912+0.0013 


1.5+0.3 




2 


4026 


2374+2383 


Fe ii* 


0.6916+0.0013 


6.0+1.2 




3 


4375 


2587 


Fe ii 


0.6915+0.0012 


0.8+0.3 




4 


4401 


2600 


Fe ii 


0.6889+0.0012 


1.1+0.3 




5 


4730 


2796 


Mg ii 


0.6913+0.0011 


2.7+0.4 




6 


4740 


2803 


Mg ii 


0.6909+0.0011 


1.2+0.3 




7 


4825 


2853 


Mgi 


0.6913+0.0011 


0.4+0.2" 






6302 


3727 


[On] 


0.6908+0.0008 




8.5+1.0 




6650 


3935 


CanK 


0.6900+0.0008 


1.2+0.3 






6711 


3969 


CanH 


0.6905+0.0008 


1.5+0.4 






7342 


4341 


H r 


0.6915+0.0007 




1.3+0.5 




8213 


4861 




0.6895+0.0006 




2.5+1.0 




8375 


4959 


[O m] 


0.6888+0.0006 




3.0+1.0 




8460 


5007 


[O m] 


0.6897+0.0006 




11.0+2.0 






Weighted mean 




0.6905+0.0010 






8 


3808 


2587 


Fe ii 


0.4720+0.0012 


2.4+0.4 




9 


3834 


2600 


Fe ii 


0.4714+0.0012 


2.5+0.4 




10 


4117 


2796 


Mg ii 


0.4723+0.0011 


1.1+0.3 




11 


4126 


2803 


Mg ii 


0.4717+0.0011 


1.1+0.3 




12 


4202 


2853 


Mgi 


0.4729+0.0011 


0.4+0.3" 








Weighted mean 




0.4720+0.0011 







*Blend of two Fe n lines. The restframe centroid is roughly estimated to be at 2380 A 



This identification is only tentative; it is however included in the list as it coincides with the redshift of this absorption system 



preliminary analysis of the spectrum of galaxy '1' led to an in- 
correct result because the contamination by the host galaxy of 
the GRB, at z = 0.691, had not been properly removed. 

3.4. Polarimetry 

In order to evaluate the Stokes Q and U parameters of the OT 
emission, and thence its polarization percentage P and posi- 
tion angle 0, we applied the method described by di Serego 
Alighieri (1997). Here P and 6 are obtained by using the rela- 
tion S(4>) = P cos 2(6 - <f>), in which S depends on the ratio be- 
tween the fluxes of ordinary and extraordinary components of 
the incident beam, and (p corresponds to the prism rotation an- 
gle. The values of P and 6 are calculated by fitting the above re- 



lation to the measurements obtained in correspondence of each 
rotation angle. Moreover, using this formalism, one obtains the 
two Stokes parameters Q and U from the values of S (0°) and 
5(45°), respectively. The intrinsic nature of the OT P-band po- 
larization is supported by the relative location of the OT with 
respect to that of field stars in the U vs. Q plot (Fig. 8, upper 
panel). 

After correcting for spurious field polarization, we found 
Qot = 0.016+0.006 and U 0J = -0.009+0.006. The fit of the 
data with the relation described above (see Fig. 8, lower panel) 
yielded for the OT a linear polarization Pot = 1.5+0.4 % and 
a polarization angle Got = 172°±8°, corrected for the polariza- 
tion bias (Wardle & Kronberg 1974). This latter correction is 
introduced because Pisa positive definite quantity, and thus at 
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Fig. 4. Optical-NIR broad-band spectra of the GRB020405 af- 
terglow at various epochs between ~1 and ~60 days after the 
GRB. Different symbol styles refer to different epochs. When 
needed, the broad-band spectra were rescaled in flux for clarity 
using the factors indicated in the figure. A spectral change lo- 
cated at (2.5+0.6)xl0 14 Hz is apparent in the first two epochs. 
Assuming a broken power-law shape, optical and NIR spec- 
tral slopes are all consistent with being 1.3 and 0.65, respec- 
tively, during the first ~10 days. In the last three epochs, the 
optical slope becomes 3.5, thus significantly steeper, and a pro- 
nounced spectral curvature is apparent between the optical and 
NIR ranges 



low S/N polarization levels the distribution function of P is no 
longer normal but becomes skewed, which causes an overesti- 
mate of the real value of P (Simmons & Stewart 1985). 

The use of standard polarization equations (e.g., 
Ramaprakash 1998) leads to results which are in good 
agreement, within the uncertainties, with those obtained from 
the method described above. In our treatment, we implicitly 
assumed that the OT polarization degree and angle did not 
vary across the four exposures based on the fact that the total 
flux varies only by 0.025 mag in this time interval. 

The spurious polarization contribution due to the Galactic 
ISM and produced by dust grains located along the line of 
sight (preferentially aligned in one direction) depends on the 
E(B - V) color excess according to the empirical relation 
(Hiltner 1956; Serkowski et al. 1975) Pism(%) = 9.0 E(B-V). 
In this case E(B - V) = 0.055, as mentioned above, so the 
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Fig. 6. VLT-Me/i>a/(+FORSl) fl-band image of the 
GRB020405 field acquired on 15 April 2002. The host 
galaxy of the GRB and galaxies '1' and '2' described in the 
text are indicated. The field size is about 25"x25"; North is at 
top, East is to the left. The diagonal solid lines indicate the slit 
position and width of the VLT spectrum obtained on 7 April 
2002 
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i5007 emission lines, at a redshift z = 0.472, are indicated 



total ISM-induced polarization is 0.5%. The fraction of this 
due to the close Galactic ISM has been taken into account and 
corrected by normalizing the OT polarization measurement to 
several field stars, as described in Sect. 2.3. Indeed, we find 
that the overall spurious polarization of the field stars is P sp = 
0.5+0.1 %. This means that it is entirely due to the Galactic 
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Fig. 8. (Upper panel): positions of field stars and of the OT 
(marked with a filled dot) in the plane of the U and Q param- 
eters not corrected for spurious (instrumental plus field) polar- 
ization. The OT is clearly separated from the region occupied 
by the field stars: this indicates that it has net intrinsic polariza- 
tion. (Lower panel): cosine fit of the /?-band polarimetric data 
of April 6. The best fit yields P = 1.5+0.4 % and = 172°±8° 

ISM, which further confirms that the instrumental polarization 
is negligible. We do not expect that the far (Galactic Halo) ISM 
contributes significantly to the OT polarization. The contribu- 
tion of the host galaxy to the total (OT+host) light is just few 
percent; therefore we do not expect that the host emission con- 
taminated significantly our polarization measurement (see also 
Covino et al. 2003). Concerning possible additional local dust 
absorption within the host (see the caveats of Covino et al. 2003 
about this issue), we will show in Sect. 4.2 that we do not find 
evidence for this, based on the analysis of the emission of the 
GRB020405 afterglow. Therefore we tend to conclude that the 
measured OT polarization is intrinsic. 

4. Discussion 

4.1. Light curves: SN vs. shock re-energization 

The optical light curves of the GRB020405 afterglow are con- 
sistent with a single power-law decay of temporal index a = 
1.54+0.06 between 1 and 10 days after the prompt event, with 



no temporal break. In the fireball scenario, this points to a 
nearly spherically symmetric expansion (Sari et al. 1998), as 
also suggested by our fits to the broad-band spectra of the after- 
glow (see Sect. 4.2). Our temporal slope is different than those 
obtained by Bersier et al. (2003) and Price et al. (2003), likely 
due to the longer time baseline of our sampling and to the lower 
host galaxy contribution resulting from our analysis. For the 
NIR light curves we find a slightly shallower temporal decay, 
a = 1.28 + 0.03. 

The red bump detected ~20 days after the GRB has been 
modeled by Price et al. (2003) with an emerging SN which 
is 0.5 mag dimmer than the Type Ic SN1998bw (z = 0.0085) 
- considered the "hypernova" prototype (Paczyhski 1998, 
Galama et al. 1998; Patat et al. 2001) - exploded simultane- 
ously with the GRB at its redshift. Similarly, Dado et al. (2002) 
modeled the light curve with a SN akin to 1998bw, but with the 
addition of intrinsic extinction. We have also applied a com- 
posite power-law-plus-SN fit to our VRIJ light curves, which 
extend over a longer time interval than that considered by the 
above authors. We first tried a fit with SN1998bw as a template, 
redshifted to z — 0.691, and then we used the light curves of 
another Type Ic SN, 2002ap (Nakano et al. 2002; Gal- Yam et 
al. 2002b; Pandey et al. 2003), located in M74 (8 Mpc) which 
has also exhibited some hypernova characteristics (Mazzali et 
al. 2002; see however Berger et al. 2002). SN2002ap showed 
a faster time evolution after maximum with respect to that of 
SN1998bw, and therefore it could possibly better describe the 
observed post-flattening OT decay. To construct light curves at 
the effective wavelengths of the GRB020405 afterglow obser- 
vations we interpolated the light curves of the two model SNe 
or extrapolated to the rest-frame UV range. In order to opti- 
mize the fit, the SN1998bw and SN2002ap temporal profiles 
were dimmed by 0.6 mag and brightened by 1.3 mag, respec- 
tively. 

In Fig. 9 are reported the power-law+SN fits for the case in 
which the SN2002ap was used as a template. Fig. 10 shows the 
fit with a power law plus SN1998bw in the R band only. Either 
using SN1998bw or SN2002ap as a template, the fit is not sat- 
isfactory in the V band. This can be expected, since this wave- 
length range corresponds to the near-UV at the GRB redshift, 
where SNe suffer severe intrinsic extinction in the envelope. 
In the other bands, the fit is marginally satisfactory for both SN 
templates. Allowing for a positive or negative time lag between 
the SN and GRB explosions does not produce any significant 
improvement in the fit. 

Alternatively, the afterglow temporal variability may be re- 
produced by both density and energy variations (e.g., Nakar et 
al. 2003). Specifically, the observed bump may be interpreted 
within the scenario of shell collisions proposed by Kumar & 
Piran (2000), in which a slower shell runs into a faster but de- 
celerating shell and re-energizes it. The amplitude of the bump 
depends on the ratio of the energy amounts carried by the two 
shells. The optical light curves of GRB020405 exhibit a larger 
temporal slope (a' = 1.85+0.15) after the bump than before 
it (a = 1.54+0.06), which would imply also a change in the 
electron energy distribution shape, an effect which is not pre- 
dicted by Kumar & Piran (2000). While this circumstance may 
limit the applicability of their picture to our observations, a 
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Fig. 9. V (upper left), R (upper right), / (lower left) and J (lower right) light curves of the GRB020405 afterglow; to indicates 
the time of the GRB onset. The data are corrected for the host galaxy contribution and for Galactic absorption, and fitted with a 
power law with a = 1.54 (dashed line) plus a SN2002ap at z = 0.691 (dot-dashed line) brightened by 1.3 mag with respect to the 
original one. The solid line represents the best-fit model. In these fits the SN was considered as occurred simultaneously with the 
GRB. Light curves of the SN2002ap template end at ~80 days after the explosion because after that date they are poorly sampled 



similar but more refined model could be consistent with them. 
Our light curves could be also qualitatively reproduced by an 
analogous scenario proposed by Ramirez-Ruiz et al. (2001), 
in which light variations are caused by the interaction of the 
blast wave with the density bumps in a pre-ejected stellar wind. 
However, this interpretation would not be self-consistent, since 
the hydrodynamics of the GRB020405 afterglow points to a 
homogeneous, rather than to a windy, circumburst medium (see 
below). 

Other views are offered by the model developed by 
Beloborodov (2003), in which the fireball is interacting with 
a trailing neutron shell, or by that in which a SN remnant, 
located around the GRB progenitor and excited by the GRB 
itself, cools down on time scales of weeks (Dermer 2002). 
Unfortunately, no thorough quantitative descriptions of the ef- 
fects induced by these interactions on the spectral and temporal 
behaviour of OTs are available at the moment. 

The dust echo model (Esin & Blandford 2000; Reichart 
2001) seems instead to be a less viable interpretation for the 



emission excess because it predicts bluer spectra at the time 
of the bump compared with those displayed by the early-time 
observations assuming scattered light dust echo; this dust echo 
emission is ruled out by the data. This model however, in case 
of thermal dust echo, also foresees a possible spectral peak in 
the observer's frame far-IR, which we cannot test. However, 
the increase in the radio flux noticed by Berger et al. (2003) 
around day 20 after the GRB may be consistent with this ex- 
planation. The microlensing interpretation (e.g., Loeb & Perna 
1998) of the phenomenon is instead unlikely, as any variation 
is expected to be achromatic in this case. 



The temporal decay steepening seen in the optical at day 
~20 after the GRB (a' = 1.85) might also be interpreted as 
a transition to a non-relativistic regime (Dai & Lu 1999) for 
the same electron distribution shape derived below in Sect. 4.2. 
However, in this case we would not expect significant spectral 
steepening after the transition, as opposed to what is observed. 
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4.2. The broad-band spectrum 

Mirabal et al. (2002) propose a wind model for the GRB020405 
afterglow expansion: the blast wave propagates in a medium 
enriched before the GRB explosion by a wind ejected from 
a Wolf-Rayet star. Under the hypothesis that the cooling fre- 
quency is higher than the X-ray energies at the epoch of their 
Chandra observations, those authors find, based only on the X- 
ray data, an electron distribution slope of p = 2.7. However, 
the X-ray spectral shape, with a spectral index j6 x = 0.7+0.2, 
is significantly flatter than the one we derive in the optical, j6 opt 
= 1.3±0.2. If one allows for some intrinsic extinction within 
the GRB host galaxy, one may recover the agreement between 
the optical and X-ray spectral indices; however, in this case it 
would be difficult to explain why the temporal decline rate is 
different at the NIR, optical and X-ray frequencies (Chevalier 
& Li 2000). If instead no intrinsic extinction is invoked, the 
flatness of the X-ray spectrum makes it impossible to recon- 
cile it with synchrotron emission, predicted by the standard 
fireball scenario to be produced as a consequence of the inter- 
action between the external shock and the interstellar medium 
(Sari et al. 1998) or the circumburst medium enriched by a pre- 
ejected stellar wind (Chevalier & Li 2000). Therefore, we fa- 
vor an interpretation of the X-rays as Inverse Compton emis- 
sion, as also proposed by Mirabal et al. (2002), and as in the 
case of GRB000926 (Harrison et al. 2001). Since the X-ray 
spectral slope is consistent, within the errors, with the NIR 
slope, the same population of relativistic electrons may have 
produced the NIR photons via synchrotron radiation and the 
X-rays via Inverse Compton scattering off afterglow radiation 
(synchrotron self-Compton). 

The optical and NIR data are inconsistent with a wind af- 
terglow model (Chevalier & Li 2000): the spectral break be- 
tween the NIR and optical bands (Fig. 4) could not correspond 
to the synchrotron injection break (v m ), since then the observed 
NIR slope would be expected to be much flatter than observed; 
therefore, this break must be due to the synchrotron cooling 
frequency (v c ). However, in this case the NIR temporal decline 
should be faster than seen in optical, contrary to what is ob- 
served. Assuming instead an isotropic and adiabatic expansion 
in an external medium of constant density ("homogeneous" 
model; Sari et al. 1998) and a slow cooling regime, the spec- 
tral slopes bluewards and redwards of the cooling frequency 
v c - 2.5 x 10 14 Hz imply an electron energy distribution index 
p = 2.5±0.3. This is consistent with the (better constrained) 
value, p = 2.7+0.1, obtained from the NIR and optical decay 
rates under the assumptions of the homogeneous model with 
no need for additional local absorption in the host. 

From this, and from the above considerations on the broad- 
band X-ray-optical-NIR spectral flux distribution, we conclude 
that the fireball emission is affected by negligible wavelength- 
dependent extinction associated with the host galaxy. Similarly, 
the OT detection in the U band (/l res tframe ~ 2100 A) suggests, 
in a model-independent way, negligible intrinsic absorption. 

We note that a collimated fireball before the jet break (Sari 
et al. 1999; Rhoads 1999) and an isotropic fireball (Sari et al. 
1998) are described in equivalent ways by the standard model 
and that they cannot be discerned by the temporal and spectral 
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Fig. 10. Same as the upper right panel of Fig. 9, but with 
SN1998bw as a template, dimmed by 0.6 mag 

indices. If the bump observed in the light curves at day ~20 
is due to an emission component independent from the after- 
glow (i.e., a SN), our data cannot rule out a steepening of the 
afterglow light curve, related to a collimation break, occurring 
more than ~10 days after the GRB. According to Sari et al. 
(1999), this would imply a lower limit to the jet opening angle 
of©>14°. 

The conclusion that the GRB020405 fireball expanded in a 
homogeneous medium was independently reached by Berger et 
al. (2003) from observations of the radio afterglow. However, 
these authors also suggest the presence of a jet break occurred 
=sl day after the GRB. While our optical-NIR light curve cov- 
erage does not allow us to confirm their finding, this interpre- 
tation would imply that we monitored the counterpart decay 
during the post-break phase of a collimated fireball. Following 
Sari et al. (1999), this would require that p = a ~ 1.5, which is 
inconsistent with our spectral results. 

4.3. Optical polarization 

Our /?-band polarization measurement is in line with the de- 
tections and upper limits so far obtained for GRB optical af- 
terglows (Hjorth et al. 1999; Wijers et al. 1999; Covino et al. 
1999; Rol et al. 2000; Bjornsson et al. 2002; Covino et al. 
2002c; Covino et al. 2002d; Covino et al. 2002e; Rol et al. 
2003). 

Jet geometries are thought to produce observable polariza- 
tion (Sari 1999; Ghisellini & Lazzati 1999), however, no indi- 
cation of a jet structure is derived from the light curve analysis 
for the GRB020405 afterglow (see Sects. 3.1 and 4.1). These 
characteristics, i.e., the presence of detectable polarization and 
the absence of a light curve steepening (at least in the first 10 
days), are reminiscent of GRB990712 optical afterglow (Rol 
et al. 2000; Sahu et al. 2000), although a jet break has been 
suggested by Bjornsson et al. (2001) for the latter. 

Our values of P and 9 are consistent with those measured 
for the present afterglow by Covino et al. (2003) in the V 
band ~2 and ~3 days after the GRB. However, our polarization 
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percentage measurement is at variance with that obtained by 
Bersier et al. (2003) from V-band observations acquired nearly 
simultaneously with ours. This cannot be explained by the con- 
tribution from the underlying host galaxy, because (as we re- 
marked in Sect. 3) this is negligible at the time of our polari- 
metric observation. Also, a strong dependence of the polariza- 
tion value on the wavelength is unlikely (Rybicki & Lightman 
1979) and, indeed, recent optical spectropolarimetric observa- 
tions of the afterglows of GRB020813 (Barth et al. 2003) and 
GRB021004 (Wang et al. 2003) show no substantial variations 
of P across the continuum redwards of 4000 A. 

For a spherically symmetric fireball, the net polarization 
should vanish if the magnetic field is ordered exactly perpen- 
dicular or parallel with respect to the expanding shock front 
(e.g., Waxman 1997). Therefore, rapidly variable effects re- 
moving the symmetry should be considered in order to explain 
the different polarization values measured by us and by Bersier 
et al. (2003), like turbulence or microlensing (Medvedev & 
Loeb 1999). In these cases substantial variations of 9 would 
also be expected (see Rol et al. 2000 and references therein), 
but they are not observed. Therefore, if real, this fast variation 
of P lacks a satisfactory interpretation. 

Thus, as also pointed out by Rol et al. (2000) for 
GRB990712, no current polarization model adequately ex- 
plains the entire optical polarimetry data set for the 
GRB020405 afterglow, that is, the strong variability of P and 
the simultaneous constancy of 8 when no break is seen in 
the optical light curve during the polarimetric monitoring (see 
however Bjornsson & Lindfors 2000). We independently an- 
alyzed the polarimetric data of Bersier et al. (2003) and, by 
applying our PSF-fitting photometric procedure, which differs 
from that based on the ISIS2 package (Alard 2000) used by 
Bersier et al. (2003), we obtain a 3-<x upper limit on the V- 
band polarization of 36%. This value is higher than, but still 
consistent with the measurement and upper limit obtained by 
the above authors. 

4.4. The host galaxy and the intervening absorption 

The ground-based and HST imaging reveal that the 
GRB020405 OT is located in a complex field (see Figs. 3 and 
6). Emission and absorption lines detected on the optical spec- 
tra of the OT indicate a redshift of z — 0.691, confirming our 
preliminary analysis (Masetti et al. 2002a) and the indepen- 
dent spectroscopic results by Price et al. (2003). An absorption 
system at the lower redshift of z = 0.472 is also detected in 
the OT spectrum acquired on April 7 (see Fig. 5), and we at- 
tribute it to galaxy '1' (Fig. 6), for which we detect [O n] and 
[O m] in emission at the same redshift in the VLT spectrum. 
At this redshift the angular separation (~2") between galaxy 
T and the OT corresponds to ~13 proper kpc, consistent with 
the halo of galaxy ' 1 ' being responsible for the absorption sys- 
tem at z = 0.472 detected in the OT spectrum. For galaxy '2', 
located at ~6" from the OT (Fig. 6), we also measured a red- 
shift z = 0.472, therefore the two galaxies probably belong to 
a foreground association or cluster. Galaxy '2' may also be re- 
sponsible for the absorption system, although its larger angular 



distance from the OT makes it a less probable candidate than 
galaxy ' 1' for the absorber. The presence of intervening absorp- 
tion systems is not uncommon in the OT spectra (e.g., Metzger 
et al. 1997; Masetti et al. 2001; Barth et al. 2003; M0ller et al. 
2002). However, the present case is the first for which the in- 
tervening absorber has been identified with imaging and spec- 
troscopy. 

The flux of the [O n] ,13727 can be used to determine 
the star formation rate (SFR) of the GRB020405 host galaxy. 
Applying Eq. 3 of Kennicutt (1998), we determine a SFR of ~3 
M yr _1 , consistent with the findings of Price et al. (2003). As 
noted by these authors, this is actually a lower limit to the SFR 
as no correction for the unknown local absorption has been de- 
termined. 

The detection of [O n], [O in] and Hp emission lines in 
the spectrum also allowed us to infer the oxygen abundance 
in the GRB host. By computing the parameter ^23, defined 
by Kobulnicky et al. (1999) as the ratio between the sum of 
the oxygen forbidden line fluxes and the Hp line flux, we find 
R23 = 6.5+2 which implies, from Fig. 8 of Kobulnicky et al. 
(1999), that 12 + log (O/H) = 8.0+0.25. This suggests that the 
GRB020405 host has an oxygen abundance lower (around one 
eighth) than solar; also, this metallicity yield is similar to those 
found for other GRB hosts (e.g., Price et al. 2002c). 

5. Conclusions 

The most important results of our ground-based observing 
campaign on the GRB020405 afterglow can be summarized as 
follows: 

- we have reported the first detection of the NIR afterglow of 
this GRB: the OT is clearly detected in JHK S bands; 

- thanks to our accurate method of host subtraction we find 
that, between 1 and 10 days after the GRB, the decay of the 
afterglow light curves is consistent with a single power law. 
Therefore we model the afterglow emission as synchrotron 
radiation in a fireball expanding in a homogeneous medium 
and within a jet whose opening angle is not smaller than 
14°. At X-ray frequencies, the contribution of synchrotron 
self-Compton radiation may dominate; 

- the addition of late-time (20-150 days after the GRB) HST 
points and VLT /-band observations indicates the pres- 
ence of a red bump in the VRIJ light curves around 20 
days after the GRB. This bump can be modeled with a SN. 
Alternatively, HST and late-time /-band data can be mod- 
eled by using a power law with decay index steeper than 
that of the early decline phase (a' = 1.85+0.15). This can 
be explained as due to a late shell collision in the fireball; 

- a /?-band polarimetry measure shows that the afterglow is 
polarized, with P = 1.5+0.4 % and 8 = 172°+8°. No quanti- 
tative model can explain the entire polarization data set for 
this OT; 

- the spectroscopic detection of emission lines indicates that 
the GRB is located at redshift z = 0.691. Two metallic ab- 
sorption systems, at z = 0.691 and at z = 0.472, are also 
detected in the afterglow optical spectrum. This latter red- 
shift coincides with that of a galaxy complex angularly 



16 



Masetti et al.: Optical-NIR monitoring of the GRB020405 afterglow 



close (2") to the host of GRB020405, indicating that the 
absorption system most likely originated in this complex. 
For the first time, the galaxy responsible for an intervening 
absorption line system in the spectrum of an OT is directly 
detected; 

- the host of GRB020405 appears to be a galaxy with mod- 
erate SFR (~3 M yr _1 ) and with subsolar metallicity; 

- while monitoring the late phases of the afterglow decay, 
HST revealed the presence of a further variable object lo- 
cated ~3" from the OT, in the southern outskirts of galaxy 
' 1' . While its nature is still to be investigated, its presence 
underscores the need for a careful validation of the variable 
candidates detected in GRB error boxes as counterparts of 
the high energy events. 

Our results point to the need of further systematic 
NIR/optical campaigns of GRB afterglows, coordinated with 
higher energy monitorings, to better constrain the emission 
mechanisms and the environmental characteristics of GRBs 
and of their lower energy counterparts. 
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